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Abstract 

In a bottom-up approach we investigate lepton-flavour violating processes r —> that are 
mediated by New Physics encoded in effective-theory operators of dimension six. While 
the opportunity to scrutinize the underlying operator structure has been investigated 
before, we explore the benefits of utilising the polarization direction of the initial r 
lepton and the angular distribution of the decay. Given the rarity of these events (if 
observed at all), we focus on integrated observables rather than spectra, such as partial 
rates and asymmetries. In an effort to estimate the number of events required to extract 
the coupling coefficients to the effective operators we perform a phenomenological study 
with virtual experiments. 



1 Introduction 


Within the Standard Model (SM) of particle physics, lepton flavour is conserved as long 
as the neutrino masses are exactly zero. The discovery of massive neutrinos and neutrino 
oscillations, however, shows that lepton flavour violation (LFV) is in principle allowed. For 
example the process r~ —* n ~7 can occur at the one-loop level if a tau-neutrino oscillates into 
a muon-neutrino within the loop that is supplemented with a charged vector boson. If this 
mechanism were the only source of LFV then the branching fractions of the r~ —>- fi~£ + £~ and 
t~ —> e~£ + £~ decay channels are non-zero, but tiny - 0 (lO -45±5 ) or thereabouts - and clearly 
unobservable. Many New Physics (NP) scenarios, however, predict much higher branching 
fractions, very roughly (9(1CF 10 ), which are at the edge of observability (see e.g. the review 
in [I]). Here and below £ denotes either a muon or an electron. 

If such a decay were to be discovered (see e.g. 1 ). as is not unrealistic given the hints on 
deviations from the SM in the lepton sector with the recent measurement^] of Rk = B(B —> 


/L/i + /i )/B{B —>■ Ke + e ) by the LHCb col 
fir) = (0.89^o; 37)% by the CMS experiment 


aboration [3| and the LFV Higgs decay B(h -+ 
4], it will be highly interesting to investigate the 
underlying interaction structure in order to disentangle possible NP models. In this paper we 
employ a bottom-up approach and treat the SM as an effective field theory, i.e. we consider 
higher-dimensional operators that consist only of SM fields and respect SM symmetries. Our 
aim is to gain qualitative and quantitative information on the couplings of these operators. It 
is obvious that such an approach will allow us to only gain direct insight into the couplings 
at the low scale - set by the tau-lepton mass - and not at the high scale of the responsible 
NP mechanisms. The implementation of the renormalization-group running of the couplings, 
which has been derived at one-loop level in (5j-7|, is beyond the scope of the current paper. 

As some of us have shown in a previous publication [ 8 ], the r —> 3£ decays are mediated 
by a handful of dimension -6 operators with possibly complex coefficients. In that publication 
a Dalitz-plot analysis was entertained in order to differentiate between radiative and leptonic 
operators. Clearly, a reconstruction of a Dalitz distribution requires a large data sample 
which obviously is hard to obtain for a very rare, lepton number violating decay. Thus, in 
order to obtain information on the type of operator that mediates the decay and on the helicity 
structure of the interaction, appropriate observables need to be defined. Obvious candidates 
are observables such as forward-backward asymmetries, which can be measured even at very 
low statistics. Our strategy is therefore to define observables of partially integrated phase 
space, which can be measured by a simple counting experiment. 

Additional information on the structure of the interaction can be gained by studying the 
decay of polarized tau leptons [9j. Such a polarization can be realized at e + e~ colliders 10 
running close to the t + t~ threshold with polarized electrons or electrons and positrons. As 
we shall show in this paper, taking into account the spin direction of the decaying tau lepton 
allows us to obtain information on the structure of the interaction, even with quite small data 
samples. 

The idea of using polarized tau leptons has also been discussed in 


11 , however, in the 


1 The LHCb measurement Rk = 0.745^o o 82 i s 2.6a away from the SM expectation Rk = 1 + 0(10 4 ) and 
challenges the notion of lepton-flavour universality. 
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context of a specific NP model. Our focus is a general analysis in order to pin down the 
structure of the relevant interaction. 

Considering the possible LFV tau decays into leptons one may classify the six distinct 
channels as (a) r“ —> £~£~£ + with £ = /i, e, or (b) r~ —> /i _ e _ e + ,T _ —> e”/i _ p + , or (c) 
t~ —> e _ e _ /i + ,r _ —» /i~p~e + . The processes in (a) and (c) involve identical particles in the 
final state, whereas (b) does not. Radiative operators contribute to classes (a) and (5), but 
not (c). For definiteness and simplicity we focus our attention solely on case (a), in which 
two external mass scales suffice, the mass of the tau lepton, m T , and the mass of the lighter 
lepton, 77i. 

The paper is organised as follows: In the next section we define the effective Hamiltonian 
with open coefficients £*, where i counts through the various operators of mass-dimension 
6 . The setup of our calculation is described in detail, including the polarisation vector of 
the initial tau and the angles characterizing the position of the polarisation vector relative 
to the decay plane. As a result the totally differential decay rate is decomposed in terms of 
trigonometric functions of these angles. In chapter [3] we integrate the differential decay rate 
over the relevant parts of the phase space. The results are rather bulky in print, so we have 
diverted them to the appendix for easier reading. Chapter [4] consists of a phenomenological 
study of the two decay channels r —> 3p and r —>■ 3e. The main question we try to answer 
is how well one could determine the Wilson coefficients of the dimension-6 operators under 
the hypothesis of having discovered a low number of events experimentally. We summarize in 
chapter [5j 


2 Calculational framework 


2.1 Operator basis 


In the following, we adopt the conventions and notation of [8]j. (For more details we also 
refer the reader to that reference.). Starting point is the most general set of dimension-six 
operators respecting the SM gauge symmetries (see 12,13)). After integrating out the weak 


gauge bosons and the Higgs field after electroweak symmetry breaking, there remain four 
purely leptonic 4-fermion operators of dimension six. Ordering the operators by the chirality 
of the involved lepton fields, they can be expressed as 


rr(LL)(LL) 


tt(LL)(RR) 

■°lept 


(LL){LL) _ _ ( RR)(RR ) 

V 2- {.^L1^T L )(£ L ^£l) , #lepf = - V ^2- (^RI/J-Tr) (V-RI^P-r) , 


( LL)(RR) 

— (hlnT L ) {£rY£r) 


{RR)(LL ) 

(i) 


where we have already singled out the r-lepton fields. In this notation, the chirality struc¬ 
ture ( LR)(RL ) can be Fierz transformed to contribute to H^ p T ^ (nRl and h[^ )(LL \ Notice 
that in order to feed terms of the form ( LR)(LR ), one would have to include dimension-eight 
operators in the SM effective field theory. Following |8|, we assume that the NP scale A is 
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sufficiently large compared to the electroweak scale, that these contributions can be neglected 
(for similar reasoning in other context, see e.g. 14 15j). In an explicit UV completion of the 
SM, the dimensionless couplings gy should be determined from a matching calculation. 


Similarly, one generates radiative operators, where, at low energies, only couplings to the 
photon held have to be considered. (Contributions from intermediate W ± or Z° bosons are 
already contained in ([Tj) .) We are then left with two more terms in the effective Hamiltonian, 


tt{LR) _ 

11 roH - 


e v 


'rad = ~^j29iri\tL<Tp,T R )F' 


/J,IS 


rr(RR) _ ^ ^ (!> ~ t- \ p/rr' 

-“rad — 87rA 2 ^ rad ^ R(J llljT L '^ 


They contribute to the r —>■ i i £ + amplitude via photon exchange, schematically 

= a=„-W“’ (h^q^TR) (h“i) . 

"---V-' ' 1 ' 

£5 m r 

(ere) ■ P) 

i&eriT 


Here, for simplicity, we used the same notation for lepton fields and on-shell spinors, q M 
denotes the momentum how through the virtual photon, and a em the usual electromagnetic 
hne-structure constant. In summary, the generic effective Hamiltonian for r —> 3£ decays can 
be written as 


TT _ tt(LL)(LL) rr(RR)(RR) jj(LL)(RR) rr(RR)(LL) MLR) MRL) 

“eff — -n lept -t- -n lept I- -n lept t- -n lept n rad n rad . 

For future convenience we shall combine the six couplings into a complex-valued vector 


( 3 ) 


e = 


(LL)(LL) (RR)(RR) ( LL)(RR ) ( RR)(LL) 

yv yv yv yv Ui 


A 2 


A 2 


A 2 


v ( LR ) a e 

>9ra.d 


A 2 ’ 77V A 2 rad ’ 77V A 2 


V {RL) 

£ad 


( 4 ) 


of mass dimension (-2). 


2.2 Spin polarization of the tau lepton 

The spin polarization of a beam of particles is typically measured in the flight direction of said 
particles. In the setup of our calculation we will define the z-axis of our lab-frame coordinate 
system to be the flight direction of the tau lepton. The reference vector for the spin orientation 
is then chosen as s M = (0, 0, 0,1) in the tau lepton’s rest frame, with s 2 = —1 and s ■ p T = 0. 
In the calculation of the squared amplitude |_M|~ we will then use the Dirac matrix 

w = ($i r + m T ) 1 + 2 75 - ( 5 ) 

to project onto tau leptons with “spin-up”. In this way the spin vector only appears 
at most linearly in the squared amplitude. Note that there are then only three linearly 
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independent invariants that one can build from and the lepton momenta in the decay 
r-(p T ) —> £~(pi) W(p 2 ) ^ + (Ps)- These can be taken as 

t = S-P2 , U = S ■ p3 , V = £ af} P “ pi pi , (6) 

and therefore the squared amplitude of the spin-up tau decay can be decomposed as 

\M t (p T ,Pi,p 2 ,P3,s)\- = gi(p r ,Pl,P2,P3) + 9t(p T ,Pl,P2,P3)t 

+9u(Pr,Pl,P2,P3)u + g v (p T ,Pi,P2,P3)v . (7) 

Similarly the decay of a “spin-down” polarized tau lepton can be calculated with the help of 
the projector ybtb = (jf T + m T )( 1 — 7 s$j/ 2 . We stress that the direction in which the spin is 
measured, i.e. the z-axis (or reference vector s M ) remains fixed. With this point of view0 the 
variables defined in (| 6 ]) remain unchanged when describing the decay of a spin-down rather 
than spin-up tau lepton. In this spin-down case one finds the same functions g± } t, u ,v from above 
that describe the squared amplitude, albeit in the combination 

\M i (p T ,Pl,P 2 ,P 3 ,s)\- = gi(Pr,Pl,P2,P3) « St(p T , Pi, P 2 , Pz)t 

-9u(Pt,Pi,P2,P3)u - g v (p T ,p 1 ,p 2 ,P3)v . ( 8 ) 

Therefore only the part 2 g\ contributes to the unpolarized decay rate. The information con¬ 
tained in the functions gt, u , v would then be lost and could not be used in the effort to unveil 
the underlying operator structure. In what follows, we will use a slightly modified version of 
the decomposition (j?]) based on angles associated with these invariants. 

2.3 Euler rotations 

The three momenta Pi,f> 2 ,P 3 of the decay product span the decay plane. In general the spin 
direction s does not lie in this decay plane, and its orientation can be described with the help 
of two angles. For example, the angle between the normal of the decay plane and s is the 
first, and the angle between the momentum p 3 of the antilepton and the projection of s into 
the decay plane the second angle. For convenience we employ the technique of Euler rotations 
instead, which also serve to define the orientation of s and the decay plane: 

We define the lab frame as the reference frame (RF) in which the tau lepton is at rest and 
the z-axis is aligned with the spin vector s = e z . We will now perform rotations until this 
plane is spanned by the new basis vectors ej and ej. Furthermore the antilepton’s momentum 
shall point directly in the e~ direction. This new RF is called the decay frame, or RF’. We 
define the following order of rotations to get from RF to RF’: 

Pi = R z (oi)R y (l3)R z (ct)px , etcetera, (9) 

Alternatively one may perform the transformation —> —s^ to find the differential decay rate of a 

“spin-down” polarized tau. 
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with mathematically positive convention: 


R z (a) 


/ cos a — sin a 0 \ 
sin a cos a 0 

0 0 1 / 


Ry{P) 


/ cos (3 0 sin /3 ^ 
0 1 0 
y — sin (3 0 cos /3 


Hence the polarization direction vector in RF’ reads 


s' = (cos cc sin /3, sin a sin /3, cos (3) . 


( 10 ) 


( 11 ) 


Note that there is no dependence on the first rotation angle 7 , which reflects the azimuthal 
symmetry of the problem. The rotations are such that the decay plane coincides with the 
x' — z' plane, so that we may parameterize 




/ b 2 sin p ^ 


/ 0 \ 

Pi = ~(P2 +P3) 

, P2 = 

0 

, P3 = 

0 



\ b 2 cos tp y 


\ b J 


( 12 ) 


where the parameter <p = 0 2 3 has a physical interpretation of the angle between p 2 and p 3 
within the decay plane, which is independent of the rotation angles a, (3 , 7 . The lengths of the 
3-vectors b 2 = | p 2 \ = \p 2 \ and similarly 63 are expressed in terms of the energies Ei as 


b 2 = 




COS 023 


El - El - El + m 2 
2 V / E 2 - m 2 \JE\ - m 2 ' 


(13) 


While the momentum of the antilepton, P 3 , is well defined, the lepton pair in the final state is 
indistinguishable if the leptons are of the same flavour, and we must be careful not to double 
count. 


For the last Euler rotation, it is sufficient to consider values of a G [0, n) to lock the decay 
plane into the x' — z' plane. Therefore, in total we restrict the values of the rotation angles as 

aG[0,7r) , 0e[O,vr) , 7 G [0, 2tt) . (14) 

Furthermore it is obvious that either p\ has a positive x' component and p 2 a negative one, 
or the other way around. The parameter 0 23 is well defined as the angle between p 2 and if 
we insist on naming that lepton’s momentum p 2 which has a positiv^jx' component, i.e. 


023 £ [0, 7r) . (15) 

(However, the squared amplitudes, which we are going to calculate, will be symmetric under 
the exchange of p\ yy p 2 , as they must for identical particles in the final state.) 

3 With this convention the invariant v = ea^sP^P^Pl 35 = £ a/37<5Pi? , 2? 3 3 s<5 is negative definite. 
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In terms of the angles and energies the invariants in (j 6 | are 


t = — b 2 sin 023 cos a sin [3 — b 2 cos 0 23 cos (3 , 
u = —63 cos f3 , 

v = —m T 6 2 6 3 sin 0 23 sin a; sin 0 . (16) 

Therefore we may decompose the squared amplitude in (J7|) in terms of the trigonometric func¬ 
tions that appear above instead of the invariants themselves. Specifically we have 

\M' 0’ = Ji(E 2 , E 3 ) + J 2 (E 2 , E 3 ) cos (3 + J 3 (E 2l E 3 ) cos a sin (3 + J±(E 2 ,E 3 ) sin a sin (3 . (17) 


This setup allows us to quite easily pick out the individual contributions from J\ to J 4 by 
folding the differential decay rate with appropriate weight functions. 


2.4 Phase space 


Besides the angles in (14) we choose two of the three energies of the final-state leptons, which 
in the rest frame of the tau satisfy E 3 + E 2 + E 3 — m T . It is a straight-forward exercise to 
show that the totally differential decay rate is then given by 


d 5 T(r3 3£) 


1 | APp 


da dcos (3 d'y dE 2 dE 3 2 m T 2567T 5 


(18) 


Since (AT | J does not depend on the angle 7 one can trivially integrate over the allowed 
range. The phase space for the energies E 2 and E 3 is nearly triangular with edges that are 
smoothened out by the light-lepton mass m. One way of expressing the phase-space boundaries 
is, for example, 


1 

2 lm 


{m T - E 2 - b 2 d 2 ) <E 3 < - (m T - E 2 + b 2 d 2 ) , m < E 2 < 


rnz 


3m 2 


2?m 


(19) 


where b 2 = \JE\ — m 2 as before, and 


m 2 — 2 m T E 2 — 3 m 2 
m1 — 2 m T E 2 + m 2 


( 20 ) 


is a function of E 2 that is near 1, except when E 2 is near the endpoint where d 2 rapidly falls 
to zero. Sometimes it is more advantageous to express these energies in Dalitz-like invariants 
Sij = (Pi + Pj) 2 , for example when addressing the momentum flowing through the virtual 
photon in the radiative operators in (J2j). In terms of s^- the energies are 


E 2 


m 2 + m 2 — Si 3 
2 m T 


E 3 


m 2 + m 2 — Si 2 
2 m T 


( 21 ) 
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Phase space 



Figure 1: Available range of the energies £ 2 and £ 3 . The shaded areas show the regions in 
which cos 023 > 0 and where cos 0 i 2 > 0 . 


For later reference we also state the phase-space region in which cos 0 2 3 > 0, which is 

1 / F , x ^ P ^ m l~ 2 m T £ 2 + m 2 
~(m T -E 2 - h d 2 ) < £3 < 2(nv _ , 


m < E2 < -(m T — m ) • 

Similarly, for cos0i 2 > 0, we have the restrictions 

m T (m T — 2 £ 9 ) + 2£^ — m 2 „ 1 


2(m r - £ 2 ) 


< £3 < -(m T - £ 2 + b 2 d 2 ) , 


( 22 ) 


m < £ 2 < - (m r — m) 


(23) 


In Fig. [l]we show the phase space for all kinematically allowed values of £ 2 and £3 mentioned 
in equation (19) as well as for the regions defined by (22) and (23). 


3 Definition of the Observables 

In this section we will use the formulae from the last chapter to define observables, which 
can be measured even with sparse data samples. All the observables will be defined from the 
coarsely sliced phase space using the angular variables, so that the observables correspond to 
partial rates and forward-backward asymmetries with respect to the angles, which also involve 
the direction of the tan-lepton polarization. 
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Table 1: Coefficients of the linear combination in (26) for a few typical weight functions. Since 

/"r>\ 

the matrices A) ’ depend on the region R in the energy phase space we leave the upper index 
c on the partial rates open. 


Partial rates combination 

W (a, cos (3) 

Cl 

c 2 

c 3 

c 4 

W F c 

Z_/ a,b ab 

1 

2 

0 

0 

1 

X)a(-^al — r^j 2 ) 

6 *(cos (3) — 9(— cos (3) 

0 

1 

0 

0 

+ r 2b — rit, — r2 fe ) 

9{l-a)~9(a~l) 

0 

0 

1 

0 

12b 2(^2 b + ^3b) 

2 9(*f-a)9{a-l) 

2 

0 

0 

v/2 


3.1 Coupling bilinears 

Since lepton-flavour violating processes are rare, we focus on observables where the available 
phase space in a, (3, E 2 , E 3 is at least partially integrated. There are many such observables, 
and a measurement of them will allow us to draw conclusions on the underlying operator 
structures. We start by considering the integration of J* over the energies E 2 and E% in some 
region “R” of the phase space. The results are bilinears in the couplings £ in Ml), and can be 
expressed in terms of 6 x 6 hermitian matrices A) ' that depend on the region R, to wit 


dE 2 dE 3 J l (E 2 ,E 3 ) = ^^A < f i) ^. 

6 


(24) 


R 


From this notation it is obvious that the doubly differential decay rate in a.,/3 is obtained by 
integration over the full energy phase space (19), 

d 2 r(A ->• u) 


m; 


da dcos (3 


6 ■ 256t r 4 


e 


A {full} + A {full} 


cos /3 + Af 11} cos a sin (3 + Ajj 1111 ^ sin a sin /3 


{full} 


C 

(25) 

from which asymmetries using a, (3 are readily calculable. For example, the difference in 
the partial rates r(cos/3 > 0) — r(cos/3 < 0) will only involve A^ 11 * 1 ^. We may express 
this difference as the convolution integral over the differential decay rate (25) with the weight 
function [0(cos (3) — 9(— cos /?)], where 9 is the Heaviside step function. In order to demonstrate 
that one can pick out each individual matrix using partial rates we note that 


r d 2 r 772 

da I dcos (3 -——— W(a, cos (3) = 


-l 


da dcos {3 


6 ■ 256tt 3 


e 


E^a 


{full} 


i— 1 


f, 


(26) 


and list the coefficients c* for a few typical weight functions in Table [TJ (Note the factor of 7r 
that has been absorbed into the prefactor.) We define the following partial rates Y c ab , which 
will form the basis of our simulated counting experiments in Section |4j obtained from disjoint 
















regions of the phase space in ck, /3, (E 2 ,E 3 ). The assignments are 


a 

( y. G 

b 

0 e 

c 

R 

1 

2 

3 

4 

[0. f) 

r 7T 7T \ 

U’ 2 / 

[7r 3n\ 

12 ’ 4 / 

[?,*) 

1 

2 

[O.f) 

11.0 

1 

2 

3 

{cos 023 
{COS 023 < 0 A 

{COS 012 

> 0 } 

COS 012 < 0 } 
> 0 } 


(27) 


Splitting the range of a into more than two regions is necessary for isolating A^\ Table [l] 
also states the linear combination of these rates that correspond to the given binning func¬ 
tions. It is clear that any observable constructed from the doubly differential decay rate ( |25| ) 
corresponds to a matrix that is such a linear combination of the four Aj fu11 ^ matrices. In fact 
only two of these structures contain information, since Ag fu11 ^ = A^' 1 "’’ = 0 due to the fact 
that the two final-state leptons are identical particles and that J 3 and J 4 are antisymmetric 
in the interchange of E\ -H- Eo. In other words, any weighted integral over (25) is just a lin¬ 
ear combination of the total decay rate and the above-mentioned asymmetry. Note also that 
Affun} = 0 is a necessary condition for our discussion on the unpolarized decay rate around 
equation ([7]) as is evident from the first line in Table [I] 


If, however, we consider only the part of the phase space (22) in which the angle between 
P 2 and p 3 is between 0 and n/2 (denoted by R= {cos0 2 3 > 0}), all four matrices Al cos ^ 23>0 ' > 
are non-zero and contribute new information. (Note that R= {cos023 < 0}, which is the 
complement to the full phase space, would then not yield more information.) Similarly the 
angle between the two leptons carrying momenta p\ and p 2 can be utilized. The region R= 
{cos 0 i 2 > 0 } contributes two more non-zero matrices as again ^ 4 h Ob@12>0 ) _ ^(cos</>i 2 >o) _ q 
W e therefore count eight different partial rates from which to construct observables. 


Since the vector of coupling constants, £, contains six complex parameters, i.e. twelve real 
unknowns, the above partial rates do not suffice to solve the system. One way out would 
be to obtain independent information on the coefficients £5 and £ 6 from the radiative decays 
r —> /zy. Alternatively, one could, of course, divide the phase space into more (i.e. smaller) 
regions R , but this would only make sense if sufficient signal events had been measured. For 
the time being, we will restrict ourselves to a handful of benchmark scenarios that will be 
defined and analyzed in Sec. [|J 

3.2 Calculation of the matrices A ( f v) 

The tree-level calculation of the squared amplitude is straightforward, and we have collected 
the various parts in Appendix [Aj The mass m of the light leptons in the final state has been 
kept finite in our calculations. The task is now to integrate the resulting functions Ji(E 2 , E 3 ) 
over (part of) the phase space R to arrive at the matrices A- R) . Their entries are functions of 
the mass ratio ratio e = m/m T , which is small. It is tempting to state the results analytically 
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in an expansion of e, which can be done using the method of regions 16 . We give a brief 


discussion of this strategy on two illustrative examples which can be found in Appendix [Bj 
For practical purposes, it suffices to evaluate the matrices numerically, and we list the results 
in Appendix [C] for both r —» 3/i and r —> 3e decays. 


3.3 Comparison with the literature 


The proposal of utilising the polarization vector of the tau for asymmetries 9] was preceded 
by an analysis within the context of the littlest Higgs model with T parity |TT] , in which 
the authors also considered lepton-flavour violating decays of polarized (anti-) taus into light 
leptons, among other channels. Contrary to our assumption stated after equation (jl| that 
(. LR)(LR ) and ( RL)(RL ) operators are suppressed by (v 2 / A 2 ), Goto et al. included them 
in their work. Furthermore only the leading-power expressions in (■ m/m T ) are considered - 
massless final-state leptons, in other words. With this approximation there is no interference 
between operators of different chirality, which require a mass insertion. In this limit the cou¬ 
plings of the above-mentioned extra operators, called g l Ls and g l Rs enter simply by adding them 
to gy L ' >( ' LL) and g {Rr< ^ RT ^ (i n our notation) in quadrature, see (A8c) of 


11 


When neglecting final-state lepton masses one encounters unregulated collinear and soft 
divergences in the contribution from the radiative operators at the boundaries of phase space, 
where the intermediate photon propagator can become on-shell. Goto et al. introduced a 
cutoff parameter 5 > 0 in order to stay away from this boundary, so that their expressions 
become functions of 5 rather than ( m/m T ). Here, in this present work, we went further in that 
the final-state lepton masses remain finite, and therefore all interference terms are accounted 
for and all soft and collinear divergences are naturally regulated. 


We have compared some of the terms in the fully differential decay rate stated in 11 


with the corresponding leading-power approximations of our results in Appendix [A] and find 
agreement after accounting for the different setups. 


4 Phenomenology 


We start with assuming some values for the LFV couplings, on the basis of which we generate 
Atotai events and bin them into A f£ b counts congruent with the definition of the partial rates 
in ( [27| . Our goal is then to reconstruct the couplings from a simple and straight-forward 
least-square fit of the binning counts, i.e. fitting the bin probabilities, 


Prot&K) = Ft = 


pA' ,ull >e 


(28) 


to the fraction of events in that bin, J\f£ b /J\l total- Note that the bin probabilities are invariant 
under a simultaneous rescaling of all couplings, £ —» (u£. In order to avoid this flat direc¬ 
tion in the y 2 function we impose a condition that breaks the rescaling invariance, to wit 
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£t^{ full }£ = ]_ _ We stress that a thus fitted result only reflects the relative strength of the 
couplings to each other, and not the values of the couplings themselves. Those can be inferred 
from the total decay width of this process and the lifetime of the tau, which is not our main 
focus hereEl 


As we have seen in Section 3.1 there are eight independent matrices from which we can 


calculate observables. One of them, ^l COf = d> 2 3> °} ^ pj-g^gg the imaginary parts of the couplings 
and does not contribute if the couplings are real. One may choose many different sets of ob¬ 
servables to fit for the couplings, but here we simply use the bins from which these observables 
are calculated themselves. The least-square fit is thus performed by minimizing the function 


x 2 (£,h) = X] 


a,b,c 


- t'Auu 


" total 




+ 2 „ 




(29) 


where /i is the Lagrange multiplier used to impose the above normalization condition. We use 
the unsophisticated statistical estimator for the uncertainty of each bin count AJ\f^ b = \/JV^ b 
if Mg b 7 ^ 0 and 1 otherwise. 


In the following, as a premature study, we are assuming only real couplings, for simplicity. 
In this case the above-mentioned matrix ^ 4 l cos ^ 2i>0 l does not contribute, and there is no need 
for splitting the range of a into four distinct bins. Hence we will combine the bins with a = 1 
and a = 2 into one single bin, as well as the bins a = 3 and a = 4. We will therefore fit 12 
observables (seven of which are independent) to 6 real unknowns. 


4.1 Muonic final state 

For a first impression we have randomly chosen a scenario with real couplings: 

Scenario a) £ T = cu (1.3, 0.4, —4.0,0.0, 2.4,1.2) 

£ t Af ull} £ = 1 =► u = 0.1086. (30) 

A typical entry in £ is therefore (9(0.1) to which we may compare the errors. From the prob¬ 
abilities we generate A/totai events that are distributed in the bin counts M^ b . These serve as 
the output of our virtual experiment. Notice that in this example, all couplings have entries 
of similar size (or happen to be zero). At this point we pretend that the couplings £ are 
unknown and proceed with the fit. The outcome £ fit will in general deviate from the scenario 
input £ m , and we may form the deviation vector A£ = £ fit — £ m , not to be confused with the 
individual fit errors h£ fit . We then repeat this virtual experiment one thousand times and 
display the entries of the deviation vector A£ in histograms, which peak around zero with a 
certain width (see Fig. [2]). If a histogram shows a normal distribution then the width (2cr) 
coincides with twice the mean of the fit errors 5£ fit . In Table [ 2 ] the resulting widths of the 

4 We stress however, that a measurement of the total decay width itself does depend on the underlying 
distribution over the energy and angular phase space, and therefore - in case of only a few events - requires 
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Tabic 2: Histogram widths 2ai as the error estimator for samples containing TVtotai events, as 
well as the mean fit uncertainty for Scenario a). 


Wtotal 


(< 72 , <55?) 

(*3,5*?) 

(< 74 , an 

K5£s fit ) 

K<5£ 6 fit ) 

20 

(0.51, 0.51) 

(0.51, 0.69) 

(0.47, 0.85) 

(0.50, 0.88) 

(0.15, 0.14) 

(0.17, 0.19) 

100 

(0.29, 0.26) 

(0.36, 0.42) 

(0.31, 0.43) 

(0.36, 0.48) 

(0.09, 0.05) 

(0.14, 0.12) 

1000 

(0.09, 0.10) 

(0.18, 0.19) 

(0.18, 0.22) 

(0.17, 0.24) 

(0.02, 0.02) 

(0.03, 0.04) 

10000 

(0.03, 0.04) 

(0.05, 0.06) 

(0.12, 0.11) 

(0.09, 0.09) 

(0.01, 0.01) 

(0.01, 0.01) 


histograms (= 2<jj) are listed for a few numbers of events A/totai- We calculate a* as the stan¬ 
dard deviation of the data underlying the histogram. The reader should keep in mind that 
the input values of the couplings are £?( ~ (0.14,0.04,-0.43,0.00,0.26,0.13) which is to be 
compared to the uncertainty of the fit results. Regardless of whether we interpret or 
as the typical uncertainty of the coupling £™, the conclusion is that with a few tens of events 
only the couplings of the radiative operators, £5 and £ 6 , can be extracted in any meaningful 
way, while the couplings of the leptonic operators, £1 through £4 require hundreds - if not 
thousands - of events. This pattern is a direct consequence of the magnitude of the entries in 
the matrices, see e.g. (fe5|): The largest entries are the log-enhanced diagonal elements 
of the radiative sector, which is why the observables are quite sensitive to £5 and £ 6 (unless 
these couplings happen to be generically suppressed in a particular class of NP models under 
consideration). 


One may ask if the particular choices of parameter values in scenario a) have any influence 
on the outcome of this sensitivity study. We therefore repeat the above procedure with the 
following twist: 

Scenario b) For each of the one thousand samples the couplings are chosen at random 
and rescaled to abide by £^A| full ^£ = 1 . (31) 


We generate a random set of couplings, drawn from a finite interval, which we choose symmetric 
around zero and universal for all couplings £,, and rescale to build a sample of vectors £. 
This procedure removes the preference towards particular values for the couplings, although 
implicitly it assumes that they are still of the same order of magnitude. Notice that, as a 
consequence for the condition (31), the resulting distributions in the sampled couplings £, 
are not flat. In Figure [3] we show the distributions for £4 and £5, exemplary for couplings to 
leptonic and radiative operators, respectively. Note that the means for the absolute values are 
roughly around 0.2 in both cases, in accordance with our previous observation that couplings 
are of order 0(0.1). We present the corresponding fit results in Table [3j which may be 


the information about the relative size of the individual couplings in an essential way. This also affects the 
experimental procedure to generate bounds on T(r —> 3^). 
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Figure 2: Example histograms of the difference A£i = (left) and A£ 2 (fight) for 

A/totai = 100 events in scenario a). The standard deviation of the left distribution is <j\ = 0.288, 
while the mean of the individual fit errors is <5^ = 0.256. This example shows a shape that 
is more peaked than a normal distribution. The right distribution has cr 2 = 0.362 which is 
less than 5^ = 0.421, and the shape is more box-like than a normal distribution. 


Table 3: Histogram widths 2<jj as the error estimator for samples containing Atotai events, as 
well as the mean fit uncertainty for Scenario b). 


A/total 


K<5£ 2 fit ) 


f) 


K^ 6 fit ) 

20 

(0.52, 0.59) 

(0.53, 0.62) 

(0.57, 0.82) 

(0.59, 0.82) 

(0.17, 

0.16) 

(0.18, 0.17) 

100 

(0.40, 0.35) 

(0.41, 0.37) 

(0.45, 0.43) 

(0.48, 0.43) 

(0.12, 

0.08) 

(0.14, 0.09) 

1000 

(0.22, 0.18) 

(0.23, 0.18) 

(0.24, 0.18) 

(0.22, 0.18) 

(0.06, 

0.04) 

(0.06, 0.04) 

10000 

(0.09, 0.07) 

(0.11, 0.07) 

(0.09, 0.07) 

(0.10, 0.07) 

(0.02, 

0.01) 

(0.02, 0.01) 


juxtaposed to the previous scenario. In general the widths a\ have increased due to the 
fact that one is convoluting the “fixed-coupling results” with the “coupling distributions”, 
leading to more pronounced shoulders in the resulting distributions. However, the mean 
fit uncertainties 8^' are roughly the same, with some entries larger than their fixed-coupling 
counterparts in Table [2j and some entries smaller. Our general conclusion remains unchanged. 

We stress again that even with these randomized couplings for each virtual experiment 
there is still a build-in assumption: that all couplings are of the same order of magnitude with 
a flat distribution. However, depending on the explicit New Physics model (see e.g. IT—221) the 
relative importance of radiative and 4-lepton operators can be quite different. We therefore 
also considered a scenario in which the radiative operators may be loop induced and are 
therefore accompanied by Wilson coefficients that could be much smaller compared to those 
of the leptonic operators. We may name this Scenario b’). To be concrete, here we assume 
a flat distribution for the couplings of the radiative operators that have a smaller support 
interval, by a factor of 1/(47 t) 2 . Again the absolute errors quoted in Table [3] remain the same, 
up to small variations. However, the corresponding distributions in the spirit of Figure [3] are 
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Figure 3: Distributions of randomly chosen couplings, rescaled for = 1, i.e. Sce¬ 

nario b). Shown here are (left) and £,( n (right). The remaining couplings of leptonic 
operators, through are distributed similar to and the second radiative coupling, 
similar to The mean absolute values of any coupling is around 0.2. 

such that is of 0(1) - representativ for all couplings to leptonic operators - and the ones 
to radiative operators are of 0( 1 / 100 ). In this case it is the leptonic couplings that can be 
determined with a few dozen events, while the radiative ones are elusive. If this pattern were 
to be observed, we could also expect a suppression of the r —* £7 decay channel. Ultimately 
a global fit with all relevant decay channels would be in order. 

Using Asymmetries to distinguish two Scenarios. In this subsection of the analysis 
the goal is to analyze the discriminating power of our approach to distinguish two different 
benchmark scenarios that may reflect the dynamical effects of some classes of NP models. Let 
us say that the dimension -6 operators ([3]) are induced by parity-violating interactions, i.e. one 
that couples only to left-handed taus or only to right-handed taus, 

Scenario c) Only left-handed taus participate 

i = u (1,0,1,0,0,1) , £Uf ull} £ = 1 , (32) 

Scenario d) Only right-handed taus participate 

£ = a; (0,1,0,1,1,0) , £Uf ull} £ = 1 . (33) 

Specifically and for simplicity we consider in scenario c) equal coupling constants to all opera¬ 
tors involving left-handed taus, while all other couplings vanish. Scenario d) is the equivalent 
setup with right-handed taus. 

There are several asymmetries one can construct from the angles a, (3, </>i 2 , </> 23 - Some of 
them have to be combined, for example the asymmetry in cos a: according to Table [l] it is 
proportional to A ( ,J V> , but for the full phase space AJ/ u11 ^ = 0, so one needs to further cut on 
the full phase space to gain information. Let us instead look at the asymmetry of the angle 
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Table 4: Central values and standard deviations for the asymmetry Ap in both scenario c) 
and d), using 1000 virtual experiments. 


Atotal 

A^Scen. c) 

Ag(Scen. d) 

Degree of Separation 

5 

0.20 ± 0.44 

-0.21 ±0.45 

0.36 

10 

0.19 ±0.30 

-0.20 ±0.32 

0.47 

20 

0.21 ±0.22 

-0.21 ±0.21 

0.67 

50 

0.20 ±0.14 

-0.20 ±0.14 

0.85 

100 

0.20 ±0.10 

-0.20 ±0.10 

0.95 


/3, to wit 

A = r(cos/3 > 0) -T(cosfi < 0) 1 \^ {K fc _ K rc 

0 T(cOS/3 > 0) ±T(cOS/3 < 0) Motai 2 - al a2 


(34) 


Given the above two models we can readily calculate the theoretical expectations of Ap = 0.202 
(Scenario c) and Ap = —0.202 (Scenario d). Even with very few total events one can indeed 
already obtain an impression which scenario to prefer, as is shown in Table |4| Here, again, 
we have repeated the virtual experiment 1000 times to produce a distribution of results from 
which we estimate the typical error as the standard deviation of the distribution. The central 
values of the distributions fluctuate a little around the theoretical expectations due to the 
finite number of virtual experiments, but even with a rather small A/"t 0 tai the likelihood for one 
scenario over the other is significant]^] The “Degree of Separation” between the two scenarios 
is calculated from the overlap of two normal distributions, p c (A) and p c (A), with the given 
central values and standard deviations of scenario c) and d), respectively, 


Degree of Separation 


OO 

1 - I dAMm[p c (A),p d (A)\ . 


—OO 


(35) 


With this definition the Degree of Separation is zero for completely overlapping distributions 
and asymptotically approaching 1.0 for distributions that are far apart. 


4.2 Electronic final state 

When the final-state leptons are electrons the operators and Wilson coefficients are different 
and in general independent from the above-mentioned decay into muons. Since electrons are 
much lighter than muons the entries of the matrices change, and the results can be found 
in Appendix [C] Again we start our phenomenological game analogous to the muonic case 

interestingly the situation is even better for the corresponding scenarios in which all radiative operators 
are absent, £5 = = 0. The central values then shift to Ap = ±0.257 and the errors remain the same, so that 

the separation between the two datasets becomes more pronounced. 
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above with 


Scenario a) 


i T = uj (1.3, 0.4, -4.0,0.0, 2.4,1.2) 


= i 


u = 0.0503 


(36) 


Although we assume here the same relative coupling strengths as in the muonic example above, 
we remind the reader that there is no relation to the muonic case, and the goal is simply to 
gain insights into our overall ability to determine the couplings from A/totai events. However, 
the normalization condition leads to typical coupling magnitudes that are about half as large 
as in the muonic case, see (30). Comparing the results in Table [5] with those in Table [2j 


we observe that the errors on the couplings of leptonic operators are somewhat smaller for 
electronic final states, but not by half as is the case for the couplings themselves. However, 
the radiative couplings are now far easier to accurately determine. For example, with only 
A/totai = 20 events the typical fit result yields £5 = 0.121 ± 0.039, whereas £1 = 0.065 ± 0.42 
does not allow us much insight. 


Scenario b) For each of the one thousand samples the couplings are chosen at random 

and rescaled to abide by ^A| fu11 ^ = 1 . (37) 

I 11 Figure [4] we show the distribution obtained from randomizing the couplings and rescal¬ 
ing. Again the typical size of couplings is of order 0(0.1). Just as in the previous case we 
observe no significant difference from our finding with fixed couplings (compare Table [5] to 


Table 5: Histogram widths 2<jj as the error estimator for samples containing A/totai events, as 
well as the mean fit uncertainty for Scenario a). 


A/total 

K^i fit ) 

( 02 , 5^) 


(cr 4 , J£f) 


K££ 6 fit ) 

20 

(0.33, 0.42) 

(0.42, 0.57) 

(0.31, 0.74) 

(0.31, 0.86) 

(0.055, 0.039) 

(0.055, 0.067) 

100 

(0.24, 0.27) 

(0.24, 0.37) 

(0.14, 0.38) 

(0.18, 0.42) 

(0.014, 0.014) 

(0.043, 0.037) 

1000 

( 0 . 12 , 0 . 12 ) 

(0.14, 0.16) 

(0.08, 0.17) 

(0.11, 0.15) 

(0.004, 0.004) 

(0.009, 0.010) 

10000 

(0.05, 0.05) 

( 0 . 10 , 0.08) 

(0.09, 0.09) 

(0.08, 0.07) 

( 0 . 002 , 0 . 002 ) 

(0.004, 0.004) 


Table 6 : Histogram widths 2a as the error estimator for samples containing A/totai events, as 
well as the mean fit uncertainty for Scenario b). 


A/total 


K^ 2 fit ) 

K^ 3 fit ) 




20 

(0.38, 0.48) 

(0.40, 0.49) 

(0.38, 0.75) 

(0.37, 0.76) 

(0.060, 0.049) 

(0.066, 0.052) 

100 

(0.28, 0.31) 

(0.29, 0.34) 

(0.25, 0.38) 

(0.25, 0.39) 

(0.037, 0.023) 

(0.043, 0.026) 

1000 

(0.18, 0.17) 

(0.18, 0.17) 

(0.16, 0.15) 

(0.17, 0.16) 

(0.010, 0.009) 

(0.010, 0.09) 

10000 

(0.09, 0.07) 

(0.09, 0.07) 

(0.08, 0.06) 

(0.07, 0.06) 

(0.006, 0.004) 

(0.005, 0.004) 
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Figure 4: Distributions of randomly chosen couplings, rescaled for = 1, i.e. Sce¬ 

nario b). Shown here are £j n (left) and £g n (right). The remaining couplings of leptonic 
operators, through £ 4 °, are distributed similar to and the second radiative coupling, 
Q n , similar to The mean absolute values of any coupling is around 0.1. 


Table [o]) . Next, we look at the ability to distinguish two scenarios, 

Scenario c) Only left-handed taus participate 

£ = w (1,0,1,0,0,1) , £Ur i} £ = 1 , 

Scenario d) Only right-handed taus participate 

£=w(o, i,o, i,i,o) , ey “y = i 


(38) 

(39) 


The central values for the asymmetry (34) in these two scenarios is Ap = ±0.387. We 
show the outcome of 1000 simulations for each A/" to tai hi Table [7J It is notable that even very 
few events suffice for the distinction of these two scenarios. 


Table 7: Central values and standard deviations for the asymmetry Ap in both scenario c) 
and d), using 1000 virtual experiments. 


-^total 

*4^(Scen. c) 

Ap (Seen, d) 

Degree of Separation 

5 

0.40 ±0.413 

-0.39 ±0.412 

0.662 

10 

0.40 ±0.280 

-0.37 ±0.298 

0.817 

20 

0.39 ±0.214 

-0.39 ±0.204 

0.938 

50 

0.39 ±0.131 

-0.38 ±0.129 

0.997 

100 

0.39 ±0.091 

-0.39 ±0.095 

1.000 
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5 Summary 

Since the decay r —>• 3£ is very clean, already a single event of this type would immediately 
imply New Physics, since the prediction of the Standard Model (extended by including the 
neutrino masses through the Weinberg operator) is practically zero. However, once such a 
decay would be observed, the nature of the underlying interaction has to be uncovered. In 
general this would require to study decay distributions, which is impossible with only a few 
events. 

In this paper we have studied the LFV decays of tau leptons into three muons or electrons, 
including a polarization of the tau lepton. Such a polarization can be generated from running 
an e + e _ collider in the vicinity of the t + t~ threshold with polarized electrons and/or positrons. 
Experimentally this could be realized at BES III, once a polarization of the beams would 
become possible. 

We have considered a general, model-independent set-up for the interaction mediating the 
t decay, which amounts to parameterizing the effective interaction in terms of a few operators. 
Any specific model would correspond to specific values for the coupling constants in front of 
theses operators, and hence even a rough measurement of these couplings could discriminate 
between different NP models. 

However, when determining the couplings in view of very sparse data samples we are 
required to define proper observables, which we have discussed in this paper. All observables 
are of the same nature as a forward-backward asymmetry, and therefore can be measured by 
a simple counting experiment. 

On this basis we have performed a feasibility study on how precisely one could assess the 
values of individual LFV couplings, based on only a small number of total signal events. It 
turns out that for some simplified cases (i.e. assuming short-distance coefficients to be real, 
or particular chiral patterns) different NP scenarios could already be distinguished with a 
quite small number of events. In the case that LFV could be experimentally established, our 
procedure could be easily extended by refining the binning for the energy phase space and by 
including independent information on radiative r —> £7 decays. 
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Appendix 


A Operator matrix elements 

Below we list contributions to the squared amplitude |A4 1 '| 2 . For legibility the symbol f 
denoting that the tan lepton in polarized spin-up in the s direction is suppressed. We start 
with the contributions from squared lcptonic operators, which are 

IM ( S )(Li, | 2 = (^■(Pl)h'(P2)/i + (P3)|^eff LKiL) | r_ (Pr)) 

= |6| 2 32(pi • p 2 ) [( p T ■ Pa) - m T (s ■ p 3 )] , 

= I6| 2 32(pi-p 2 ) [(p T -p 3 ) +m T (s -p 3 )] , 

\M™**>\* = |6| 2 8 {(Pi • p 3 ) [( Pt ■ P 2 ) ~ m T (s ■ p 2 )] + (p 2 • Pa) [( Pt ■ Pi) ~ m T (s • pi)] 

- m 2 [(p T ■ p 3 ) - m T {s ■ p 3 )} | , 

\ M (RR)(LL,\ > = |^ 4 |2 8 |( Pi . P3 ) [( Pr . p 2 ) + m T (s • p 2 )] + (p 2 ■ p 3 ) [(p T ■ Pl ) + m r (s • pi)] 

— m 2 [(p T • p 3 ) + m T (s • p 3 )] | . (40) 


Next, the interference terms from lcptonic operators, which are 


>C (i 1 , aC i,rr, ‘ 

U (LL){LL) k ARR){LL)* 
• /Vl lept yw lept 

Lpf 1 (fiR) Lpf (RR) * 

<? ,RR| aC ,,11) ' 

ac mrr x= r ? )(ii) * 


(-32)m r m 3 , 

£i£a 8m 2 [m T (s • p 3 ) - (p T • Pa) + nr 2 ] , 
8m r m [m T (s • p 3 ) - (p T • p 3 ) + m 2 ] , 
6^3 8 m T m [~m T (s • p 3 ) - (p T • p 3 ) + m 2 ] , 
8 nr 2 [—m T (s ■ p 3 ) - (p r • p 3 ) + m 2 ] , 
6^4 &rn T rn [(p x • p 2 ) - 2m 2 ] , 


(41) 
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as well as the complex conjugate expressions. The contributions from radiative operators read 


l^ad^l 2 = I£ 5 1 2 2 m 2 j — 6m 2 


2 m 


Zj / / 6 . o cy. / o , v o\ 

1 — 2 — (m T — m ) (m r — 2m T (s • p 2 ) — w- j 


’13 


2m : 


1 IL , n o. / o , , o\ 

1— 2 — (m T — m ) (m T — 2m T {s ■ p 1 ) — m J 


’23 


-S13S23 

1 


(m 4 — 3m 2 m 2 + 2m 4 ) (m 2 + m T (s ■ pi) + m T (s • p 2 )) 


H- 

«13 L 


1 

H- 

^23 L 


| M 


(RL) 12 
rad I 


ra 4 + 3ra 4 — 2 (m^ + 3m 2 )5 2 3 + 25^3 
+4m r (2m 2 — s 23 )(s • pi) + 2m T (3m 2 — 2 m 2 ) (s • p 2 ) 
m 4 + 3 m 4 — 2(m 2 + 3m 2 )si 3 + 2s 2 3 
+4m T (2m 2 — si 3 )(s • p 2 ) + 2m r (3m 2 — 2m 2 )(s • p 3 ) 
| £ 6 1 2 2m 2 1 - 6m 2 

2 777 ^ 

H 2~ ( m r — m2 ) ( m r + 2m T (s • p 2 ) — m 2 ) 

S 13 

2 m 


\ —2~K — m2 ) ( rri T + 2m T (s ■ pi) — m 2 ) 


3 23 


-S13S23 

1 


(m 4 — 3m 2 m 2 + 2m 4 ) (m 2 — m T (s • p 3 ) — m T (s • p 2 )) 


+ — 
«13 


mUPmSP = 


m 4 + 3m 4 - 2(m 2 + 3m 2 )s 23 + 2s2 3 
—4m r (2m 2 — s 23 )(s ■ pi) — 2m T (3m 2 — 2m 2 )(s • p 2 ) 
m 4 + 3 m 4 — 2(m 2 + 3m 2 )s i3 + 2s 2 3 

— 4m r (2m 2 — si 3 )(s • p 2 ) — 2m T (3m 2 — 2m 2 )(s • pi) | , 

. _ , , 3 N f, n o { 1 1 \ m 2 (m 2 - m 2 ) 

^ 5^6 (—4m^m) 4 + 2m 2 — + —- 

l \Sl 3 S23 J <Sl 3 S 23 


1 

H- 

^23 


(42) 
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The mixed lcptonic radiative contributions are the only ones that contribute to J 4 . We use 
the convention 

tr [YYl P l a lb] = -&£^ pa • (43) 

Since this trace is purely imaginary, we will pick up a dependence on Irn(£j£*). We find 


4m r |2m T (m 2 - s 13 - s 23 ) 

2/2 2\ / 1 1 \ 

+m T m im — m ) I-(- 

\ S 13 s 23/ 

+ 2 (m 2 + m 2 - si 3 - s 23 ) ( —-— ) i^a^sP^plp'ls 6 ) 

V-S13 S 23/ 

+ (m 2 - m 2 )(m 2 + m 2 - s 2 3 )^—^ 

S13 

■ / 2 2\ / 2 , 2 \ P 2 ) 

+(m T — m )(m T + m — Si 3 )- 

-§23 

+ [(m 2 - si 3 ) 2 - 2m 2 si 3 + m 4 ] — — 

S23 

+ [(m 2 - s 23 ) 2 - 2 m 2 s 23 + m 4 ] —— 

S13 

- (m 2 + s 23 - m 2 ) (s • Pi) - (m 2 + S 13 - m 2 ) (s • p 2 )} • (44) 

Next, 

= &>£> 4m T m|2 (m 2 - s i3 - S 23 ) + 4m T [(s • pi) + (s • p 2 )] 

2 2 
777- 777 'l 

H-[m 2 — m 2 + 2m T (s • p 2 )] H-[m 2 — m 2 + 2m r (s ■ pi)] > .(45) 

>§13 >s 2 3 J 
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Next, 


Next, 


m^m^ rv 


^ 5 ^ 32 m r |m T (s 13 + s 23 - 6m 2 ) 

2 / 2 2\ ( 1 1 
+m T m (m — m )-1- 

V-Sis S 2 3 

-2- i{£ a p l5 p 2 p^ T s ) 

\-Sl3 523/ 

— (m 2 — m 2 ) —(s • pi) — (m 2 — m 2 ) — (s • p 2 ) 
$13 -S 23 

-S 2 3(S - Pi) - Si 3 (s -p 2 ) 


- (m 2 - si 3 )si 3 + m 2 (2m 2 + Si 3 ) 


(gygi) 

-S23 


- [(m 2 T - s 23 )s 2 3 + m 2 (2m 2 + s 23 )] J. 


Sl3 


A ALR) k ARR)(ll y 

JVl rad • /Vl lept 


^ 5 ^ 4 4m r m|m - s 13 - s 2 3 

m 2 . o 2 \ ( 1 1 

-^('"V - rn ) — + — 

2 V S 13 s 23 


-2m r (— -—) i^a^sP^PsPl-S 8 ) 

\ s 13 $23/ 


m T (m 2 — m 2 )- 

(« - Pi) 

— m T (m 2 — m 2 ) 

Sl3 

m T (m 2 - si 3 )- 

(« - Pi) 

- m T (m 2 - S 23 ) 


$23 


(s • p 2 ) 

S 2 3 

> • P2) 
Sl3 


(46) 


+m T [(s -pi) + (s -p 2 )] 


(47) 


Now we mi x the second radiative operator into the leptonic ones. The expressions are the 
same as we found before, except for the signs in front of the spin-vector products (s • p\) and 
(s -p 2 ). 


M [*a M \ept )[LL) = 6^4m r m|2 (m 2 - s 13 - s 23 ) - 4 m T [(5 • pi) + (s • p 2 )] 

2 2 

+ — I'm 2 — ml — 2 m T (s ■ p 2 ) 1 + — \m 2 — m 2 — 2 m T (s • pi)l \ .(48) 
-S13 s 2 3 > 
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Next, 


• M ifd i) - A/f ie^f ' {RR) * = &£4m T {2m T (m 2 T - s 13 - s 23 ) 


Next, 


2 / 2 2s I 1 1 

+m T m (m — m )-1- 

V-Sis «23 


+2 (ml + m 2 - S13 - S 23 ) ( — -— ) i^afiySP^PrS 6 ) 

V s 13 S 23 / 


(ml — m 2 )(m 2 T + m 2 — s 23 ) 


(ml — m 2 )(m 2 T + m 2 — s i3 ) 


(s -pi) 


(s ■ P2) 


[(m 2 T - si 3 ) 2 - 2m 2 si 3 + m 4 ] 
[(ml - s 23 ) 2 - 2/n 2 s 23 + m 4 ] 


(s-Pi) 


(s ■ P2) 


+ {ml + s 23 - m 2 ) (s ■ Pi) + {ml + s i3 - m 2 ) (s 


= ^4 m T m\m 2 - s 13 - s 23 


m 2 ,2 2 \ ( 1 1 \ 

- m > W + sj 

- 2 m T (— - —\ i(£a^sP2P3Pls 5 ) 

\ s l 3 s 23/ 

, (2 2\ ( S ‘ Pi) | (2 2\ ( S ’ P 2 ) 

+m T (m — m ) -b m T (m — m ) - 

s 1 3 -s 23 

+m T (m 2 T - gi 3 ) ^ S Pl - + m T (m 2 T - s 23 ) ^ S P ^ 

$23 s l 3 


— m T [(s ■ Pi) + (s ■ 
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And finally 


Ae ) AC ,(ii) ' 


+m T m 2 (m 2 


•S13 + S23 - 6 m 2 ) 



o / ^23 ^13 \ •/ a 3 1 S\ 

-2 - ^(£a3^5P 2 P3Pr S ) 

\ s 13 ^23 / 

+ (m 2 - m 2 ) —(s ■ pi) + (m 2 - m 2 ) —(s • p 2 ) 
>Sl3 S 23 

+S 2 3(s -Pi) + Si 3 (s -p 2 ) 

+ [(ml - Si 3 )si 3 + m 2 (2m 2 + Si 3 )] ^ Pl ^ 

S 2 3 

+ [(m 2 T - s 23 )s 23 + m 2 (2m 2 + s 23 )] ^ 

-§13 ) 


( 51 ) 
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B Analytic Calculation of the Matrices A- R) 

( 

In the following, we show two examples how to obtain analytic results for the matrices A\ 
as an expansion in terms of the small parameter e = m/m T , keeping terms up to 0(e 2 ). 


Integrating a constant over the full phase space. We integrate the constant 1/m 2 
over the energies E 2 and E 3 . At leading power in e there is no problem and the result is 
1/8. However, at subleading power divergences appear at the border of the phase space. We 
regulate these divergencies by manually introducing 


1 

m 2 



V l4 ) rn 2 ’ 


(52) 


and taking the simultaneous limit 771,2 —» 0 in the end. The hrst integration - over E 3 , say 
is straight forward. We then distinguish three regions. 

• Treating E 2 as 0(m T ) yields 


f hard 


m 


8 2777/771 


2 2 r 

m m 1 

2777/772 777 / 


Vi 


t I x T ~~~ ) hr — — 


777 T 


27/2 / 2/71 


singular 


Vi 


1 + — 1 In 


h2 


(53) 


• When treating E 2 ~ m one needs to integrate E 2 G [777, 00] . This gives 

11 777 772 m T " 

—I— In- 1 -In — . 

4 2 2/7! 77! n 2 _ 

singular 


m 


^soft E2 


2777 / 7 /] 


+- 


777 

777 / 


(54) 


• A 2 is near its maximum value is akin to treating S 13 ~ m 2 and integrating S 13 G 
We find 

1 777 7/1 m T ' 

—h In- 1 -In- . 

2 jl 2 2-7/2 2/7i _ 

singular 


^soft S13 


m 2 m 2 

2-777/772 -777/ 


[4-7T7 2 , (X)]. 
(55) 


The singular terms cancel in the sum of these contributions, and the dependence on the 
auxiliary scales jii drops out as well, resulting in 


1 3777 2 777 

lhard T Jsoft E2 T ^soft S 13 / 1 X y lu 

8 2 / 77 / 


- 777 t 


(56) 
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Integrating singular pieces from the radiative contributions. Much more challeng¬ 
ing contributions arise from the radiative operators, where the intermediate photon propagator 
leads to integrands proportional to inverse powers of the invariants s t ,j. The most singular ex¬ 
pression in (42), for example, is 

1 1 


+ 


S 13 S 23 (Sl3 + ^23)^13 ( s 13 + $ 23)323 

For definiteness we are therefore focussing here on the integral 

.2 r 


S(m/m T ) = 


dEo dE 3 


mz 


phase space 


(S13 + S23) s 13 


dE 2 dE 3 


mz 


phase space 


(2 m T E 3 + 2m 2 )si 3 


(57) 


(58) 


which has a soft E 3 and a collinear S 13 singular structure. We may identify four different 
regions to solve this integral. 


1. The soft E 2 region, where E 2 ~ m and b 2 can be resolved. This region only becomes 
relevant for power corrections. 

2. The collinear S 13 ~ m 2 region, where E 2 is nearly at the endpoint, namely of the order 
0{m 2 jm T ) away from its maximum. Here d 2 can be resolved. 

3. The soft E 3 region, where E 2 is of order 0(m) away from its maximum. 

4. The bulk. 


We aim to distinguish these regions in a single, one-dimensional integral over a variable that 
differentiates between the four regions in terms of its power counting. To this end we perform 
the first integration over E 3 without approximation and are left with a single integral over 
E 2 . Let us now define E 2 as the energy variable that is smallest for maximal E 2 , i.e. E 2 = 
m 2 ~ 3m - E 2 . Finally we define the variable p = E 2 /E 2 , which will accomplish our goal: After 

substituting the E 2 = mT ^ m jyy integration variable for p E [0, ^ — 1 — W e find the 
following in the above four regions and name them 

1 . p ~ m T /m, label “uhard”, 

2 . p ~ m 2 /m 2 , label “usoft”, 

3. p ~ m/m T , label “soft”, 

4. p ~ 1, label “hard”. 

Finally we need to choose a regulator common in all regions. While p v —^ 1 is quite suitable 
we find that the regulator 

/ + \Ai 3 - 4m 2 m 2 \ v / EA 271 v-^o 1 
V \/si3 - \/si 3 - 4 m 2 l4 ) \d2j 

is more advantageous as it simplifies the practical calculation. We find 


( 59 ) 













Usoft region: Expanding the integrand and integrating 0 < p < oo we find 


Sus 4 if + 2r i 


In 


m T m 2m 
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2 fj,ifj ,2 rri 2 


+ finite. 


( 60 ) 


Soft region: Expanding the integrand and integrating 0 < p < oo we find 




1 1 
2rf + 2rj 


In 




rri'irn 


+ 


2 ml 


+ finite. 


Hard region: Expanding the integrand and integrating 0 < p < oo we find 


9 = J- + 1 

h 4?7 2 2r/ 


In 


2 2 
m 

-— H- 

2 hi /i 2 nr 2 


+ finite. 


(61) 


(62) 


Uhard region: Expanding the integrand and integrating 0<p< ^ — 1 — we find 


-S'ufe = -j -—j + finite. 
4r/ rrd 


(63) 


Again all the double and single divergences cancel in the sum of the regions and we obtain a 
leading double logarithm free of any regulators, to wit 


S(m/m T ) = 


1, 9 m 

7r 2 

9 

m 

'9 

5. 

m 1 

- In 2 — 


+ —— 


-m — 

4 m T 

~ 16 

m 2 

_8 

4 

m T \ 


(64) 


Armed with the techniques discussed above we are able to state some analytic results. We 
find in terms of e = m/m T 


Af ull} = 


/1 - 24e 2 


0 


4e 2 


—2e 2 - 48e 2 -4e 

0 

1 

— 24e 2 


—2e 


4e 2 —4e 2 - 48e : 

4e 2 


—2e 

1 

2 

- 14e 2 


e 1 - 24e 2 -2e 

—2e 


4e 2 


e 

1 

2 

- 14e 2 —2e 1 - 24e : 

2 - 48e 2 


—4e 

1 

— 24e 2 


—2e a£ ull} -126 

V -4e 

2 

— 48e 2 


—2e 

1 

- 24e 2 —12e a£ uI1} 


+ 0(€ 3 


(65) 


where 


J 55 Ul1 ^ = a 66 l1 ^ = — (8 In e + 11) + [12 In 2 e + 36 In e + 50 — 37r 2 ] 


■ 21 6 2 


( 66 ) 
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C Numerical Results for the Matrices AS 

L 

For muons in the final state we use m = 105.658 MeV and m r = 1776.82 MeV, which 
corresponds to e = 0.0595. 


A\ m ~ 


0.919 -0.004 0.012 -0.104 1.838 -0.217 

-0.004 0.919 -0.104 0.012 -0.217 1.838 

0.012 -0.104 0.453 0.050 0.919 -0.109 

-0.104 0.012 0.050 0.453 -0.109 0.919 

1.838 -0.217 0.919 -0.109 11.635 -0.628 

-0.217 1.838 -0.109 0.919 -0.628 11.635 


— 


0.889 0 -0.006 -0.103 1.461 -0.210 

0 -0.889 0.103 0.006 0.210 -1.461 


-0.006 0.103 -0.159 
-0.103 0.006 0 


0 -0.848 0.185 

0.159 -0.185 0.848 


1.461 0.210 -0.848 -0.185 -4.344 

-0.210 -1.461 0.185 0.848 0 


4.344 


For the part of the phase space where cos 023 > 0 we find 

/ 0.245 -0.001 0.003 -0.019 0.572 -0.041' 


I {cos <^>23 >0} 


0.001 0.245 -0.019 0.003 -0.041 0.572 

0.003 -0.019 0.081 0.015 0.175 -0.020 

-0.019 0.003 0.015 0.081 -0.020 0.176 

0.572 -0.041 0.175 -0.020 4.320 -0.139 

-0.041 0.572 -0.020 0.175 -0.139 4.320 


I {cos 023 >0} 


0.235 0 -0.001 -0.018 0.417 -0.038 

0 -0.235 0.018 0.001 0.038 -0.417 


-0.001 0.018 0.009 

-0.018 0.001 0 


0 -0.145 0.041 

-0.009 -0.041 0.145 


0.417 0.038 -0.145 -0.041 -1.831 

-0.038 -0.417 0.041 0.145 0 


1.831 


30 




( 

0 

0 

0 

0 

-0.276 

-0.0005 \ 



0 

0 

0 

0 

0.0005 

0.276 

{COS</>23>0} ^ 


0 

0 

0.037 

0 

0.162 

-0.036 

3 — 


0 

0 

0 

-0.037 

0.036 

-0.162 



-0.276 

0.0005 

0.162 

0.036 

-1.489 

0 


V 

-0.0005 

0.276 

-0.036 

-0.162 

0 

1.489 ) 


( 

0 

0 

0 

0 

0.253i 

0 \ 



0 

0 

0 

0 

0 

0.253i 

{cos </> 23 > 0 } ^ 


0 

0 

0 

0 - 

-0.162i 

—0.034i 

■4 — 


0 

0 

0 

0 - 

-0.034i 

—0.1627 



—0.253i 

0 

0.162i 0.034f 

0 

0 


V 

0 

—0.253i 0.034i 0.1627 

0 

0 / 


and for R= {cos0i2 > 0} we find 


^{cOS012>O} 


a {COS 012>O} 

A 1 


/ 0.083 

-0.001 

0.002 

-0.032 

0.100 

-0.065 \ 

-0.001 

0.083 

-0.032 

0.002 

-0.065 
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0.002 

-0.032 
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0.002 

0.272 

-0.033 

-0.032 
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0.002 

0.121 
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0.272 
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1.001 

-0.141 

Y —0.065 
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0.272 
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1.001 j 
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0 
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-0.032 
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-0.064 \ 

0 
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0.032 

0.002 

0.064 

-0.084 
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0.032 
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0 

-0.269 

0.037 

-0.032 

0.002 

0 

0.098 

-0.037 

0.269 

0.084 

0.064 

-0.269 

-0.037 

-0.793 

0 

y —0.064 

-0.084 

0.037 

0.269 

0 
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a {COS </*12>0} 


^{cos 012 >0} _ Q 


(69) 
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Now we repeat the calculation for electrons in the final state. Here m = 510.9989 KeV. 
This means that e = 0.0002876 is much smaller and we can observe the structure of the 
matrices much easier. We find 





^{COS 023 >o} 


/l.OO 
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2.00 
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1.00 
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A (cos 023 >O} 
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/ 0.27 
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0 
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0 
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0 
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i {cos </> 23 > 0 } 


0 

0 

0 

0 

-0.63 

0 

0 

0 

0 

0 

0 
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0 

0 
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0 
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0 
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I {cos 023 > 0 } 
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I {cOS 012 >0} 


0.08 0 0 0 0.09 0 

0 0.08 0 0 0 0.09 

0 0 0.14 0 0.32 0 
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